Exciton states in cylindrical nanowires 
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The exciton ground state and excited state energies are calculated for a model system of an 
infinitely long cylindrical wire. The effective Coulomb potential between the electron and the hole 
is studied as function of the wire radius. Within the adiabatic approximation, we obtain 'exact' 
numerical results for the effective exciton potential and the lowest exciton energy levels which are 
fitted to simple analytical expressions. Furthermore, we investigated the influence of a magnetic 
field parallel to the nanowire on the effective potential and the exciton energy. 
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I. INTRODUCTION 

Nowadays, semiconductor nanowires with a small diam- 
eter can be realized in a fast and cost effective manner. 
One of the possible growth techniques, the VLS (vapor- 
liquid-solid) technique, is a bottom-up technique which 
appears to be very promising^. In this method, a liquid 
metal cluster acts as a catalyst where reactants dissolve, 
which leads to the direct growth from the supersaturated 
alloy droplet. This growth method results in a freestand- 
ing nanowire that is not surrounded by any other mate- 
rial, in contrast to the embedded wires. 
From a theoretical point of view such semiconductor 
nanowire structures are very interesting for understand- 
ing the role of dimensionality on physical properties. Fur- 
thermore, they provide a great potential for several appli- 
cations, such as transistors 2 , diodes 3 , memory elements*, 
lasers*)*, chemical and biological sensors?, etc. 
These nanowires have a large aspect ratio and can be con- 
sidered one-dimensional structures. They have a length 
and a diameter in the micrometer resp. nanometer range. 
Electrons and holes are strongly confined in two direc- 
tions and there is no confinement in the direction of the 
axis of the wire. Here we will investigate the exciton 
binding energy and its lowest excited states. Excitons 
form the basis of optical properties reflecting the intrin- 
sic nature of low-dimensional systems. As a first step we 
investigate the potential and energy of the electron-hole 
pair. Because our aim is to present numerical 'exact' re- 
sults we will limit ourselves for the moment, to the model 
system in which the difference in dielectric constant of the 
nanowire material and its surrounding is neglected. Only 
in this special case we are able to present analytical re- 
sults. Although for 'stand alone' nanowires the dielectric 
mismatch in the wire (difference in dielectric constant 
of the material of the wire and the environment) is im- 
portant, the study of it will be postponed for further 
investigations. 

In previous work*, a similar model system was studied, 
but only the ground exciton state was obtained and no 
systematic study was presented of the dependence on 
wire radius, mass of the carriers, magnetic field and of 
the excited states. Also impurity states - both with and 



without a magnetic field parallel to the wire axis - were 
considered in previous papers***, with the charged im- 
purity fixed to the cylinder axis and with averaging per- 
formed only over the electronic motion. Furthermore, the 
effective potential for the ground state configuration was 
calculated before*^ and perfectly fits our results. But 
the aim of the present paper is to give a more general 
result with different analytical expressions for the effec- 
tive Coulomb potential that are fitted to the numerical 
results. These fits haven't been given before and provide 
us with a tool which is extremely useful for further nu- 
merical calculations. Furthermore, the binding energy of 
the exciton is calculated numerically for arbitrary wire 
thicknesses, while the results of Ref. were limited to 
thin wires. Finally, the addition of a magnetic field is 
another extension which hasn't been studied before in 
detail. The effective potential was calculated before for 
embedded wires*^ for magnetic fields perpendicular and 
parallel to the wire axis. The techniques used in this 
paper are the same as in our paper, but in the paper 
the focus is not on the effective potential, it is on the 
calculation of the exciton energy and comparison with 
the experiment. For high magnetic fields the analytical 
expressions as obtained for parabolic confinement*, can 
be accurately fitted to our results. Note that our cal- 
culations are very general, since they apply to arbitrary 
materials if they can be described within the framework 
of the effective mass theory. 



This paper is organized as follows. In Sect. [D] we con- 
sider the single particle states of the electron and hole. 
In Sect. Hill we calculate the effective potential of the ex- 
citon, using the adiabatic approximation. The numerical 
results for the effective potential as function of the inter 
particle distance are fitted to simple analytical expres- 
sions. The exciton energy is calculated numerically and 
semi- analytically in Sect. lIVl The influence of a magnetic 
field will be investigated in Sect. IVl and our conclusions 
are presented in Sect. I VII 
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II. SINGLE PARTICLE STATES 



For completeness we first consider the eigenfunctions 
and corresponding eigenvalues of a single particle in a 
circular— potential well with radius R. Adopting the 
effective mass approximation we write the Schrddinger 
equation as follows: 



2m 



V 2 V>(r) + V(f)ip(r) = Eip(r), 



(1) 



where we approximate the confinement potential of the 
wire by a hard wall potential 



V(r) 



if p < R 
oo if p > R 



(2) 



where p = y 1 x 2 + y 2 . For the freestanding wire the po- 
tential height is given by the work function which is typ- 
ically a few eV and therefore very large compared to typ- 
ical confinement energies. The wire is oriented along the 
z— axis. Since the motion in the z— direction (free par- 
ticle in that direction) is decoupled from the motion in 
the (x, y)— plane and because of circular symmetry, we 
have ij>(p, 9, z) = F(p)e~ e lkzZ , where 9 is the polar an- 
gle. The Schrodinger equation is then reduced to the one 
dimensional (ID) equation 



p 2 F(p)" + pF(p)' + ((k 2 k 2 z )p 2 f)F(p) = 0, (3) 

with k 2 = 2mE/h 2 . The solution to this equation are 
the well known Bessel functions and taking hard wall 
boundary conditions into account, one gets for the single 
particle wave functions 



1>n,i*.M*) = Cnje-Mjtf^p^e^*, (4) 
and the corresponding eigenenergies 



21.2 



E,, 



K l k 



i 

2mR? 2m ' 



(5) 



where C n ,i is the normalisation constant and /3 n i is the 
n th order zero of the Bessel function Ji(x). The zeros 
of the Bessel functions are known and can be found in 
the literature^: f3 0fi = 2.4048, /3 ,i = 3.8317, /3 , 2 = 
5.1356, j9i, = 5.5201, /3i, 2 = 8.4172', fa, 2 = 11.6198,.... 
Notice that the first excited state is |0, 1, k z ) and the sec- 
ond one |0, 2, k z ). 



H 



r, 2 



2m e 
+Vc (f e - f h ) 



V e (x e ,y e 



2m, 



-V 2 h + V h (x h ,y h ) 



(6) 



where m e (nth) is the effective mass of the electron (hole), 
V e (Vh) is the confinement potential of the electron (hole) 
and 



Vc (f e ~ fh) = 
1 



r,(7) 



47re y/{x e - x h ) 2 + (y e - y h ) 2 + (z e - z h )' 

is the Coulomb interaction potential between the electron 
and the hole. In this case, it is convenient to separate 
the motion of the two particles into that of the center 
of mass of the system and the relative motion, thereby 
introducing the quantities: Z — (m e z e + m/ l z/ l )/(m e + 
mh),z = z e — Zh,M = m e + mh,p = (m e mh)/(m e + mh). 
Now the Hamiltonian can be rewritten as 



H 



-^l e , Ve + V e {x e ,y e ) 



2m e 6Tn,h 

h 2 d 2 h 2 d 2 

2MWZ ~ 2p\Wz 



+V h (x h ,y h ) - 



1 



47re y/(x e - x h ) 2 + (y e - y h ) 2 + z 2 



(8) 



In the case of strong lateral confinement, the exciton mo- 
tion along the wire is decoupled form the lateral motion 
of the particles and we may decouple the wave function 
as follows 



*(x e , y e , x h , yh, z, Z) = e 1 c/)(z)ip(x e ,y e )ip(xh,yh)- 

(9) 

Because the Coulomb energy in the considered case is 
much weaker than the single particle confinement energy, 
we adopt the adiabatic approximation, thereby taking 
ip(x, y) as the above obtained single particle states. Mul- 
tiplying Eq. © with i/j*(x e ,y e )ip*(xh,yh) from the left 
and with ip(x e ,y e )ip(xh,yh) from the right, integrating 
over the lateral coordinates we reduce the original 6D 
Schrodinger equation for the exciton problem to a ID 
effective Schrodinger equation for the relative exciton co- 
ordinate 



III. EFFECTIVE EXCITON POTENTIAL 

A. Adiabatic approximation 

The appropriate Hamiltonian for an exciton in the effec- 
tive mass approximation is given by 



E e + E h 



h 2 K 2 
2M 



| V ' +y effW 



0) = E iot (j>{z), 
(10) 

where E e (Eh) are the single-electron (hole) energies, and 
the effective exciton potential, V e g, is given by 
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FIG. 1: Effective exciton interaction potential as a function of 
the electron-hole separation z for three different states of the 
electron and hole (Full curve: n e = l e = and uh = lh = 
for the hole; dashed curve: l e = l h = 1 and n e — n,h = 0; 
dotted curve: l e = lh = 2 and n e = rih = 0). The regular 
Coulomb potential, 1/z, is also plotted (dashed-dotted) for 
comparison. Inset: The effective potential for different sets of 
quantum numbers for the electron and hole. 



FIG. 3: Plot of the numerical data and the two Pade approx- 
imations, Pq(z)/Pi(z) and P\{z)/ P%{z), for (a) the ground 
state and (b) the first and second (inset) excited state. 

where V^ff(z) = V g q(z) / Eq and z is in units of R. The 
coefficients are given by 



FIG. 2: Density for the electron and hole in the nanowire for 
the lowest three single-particle states. 



V e&& = j j jj dx e dx h dy e dy h 

x \lp7i : l,kA x e,y e )\ 2 \^nd,kA x h,yh)\ 2 ^ 
\J (x e ~ X h ) 2 + (y e - y h ) 2 + Z 2 

Introducing dimensionless units (i.e. length expressed in 
units of R) and a — a* B /2R with a* B = 47re?i 2 //ze 2 the 
effective Bohr radius, we finally get 



-aVi 



V eS (z) cf>(z) = E c <Kz), 



(12) 



where z is given in units of R, V = V/Eq, Ec — Eq/Eq, 
E a = e 2 /47rei? and E c = E tot - E e - E h - h 2 K 2 /2M. 
To study the bound energy states of the exciton, we 
first have to determine the exciton potential. Examining 
Eq. (|1I|I . it is clear that evaluating this integral numeri- 
cally is not straight-forward because of the I jr singular- 
ity for r — > 0. 

The method we used to get rid of this problem, is ex- 
plained in Appendix El The result of the calculations 
can be found in Fig. ^ where we show the effective 
Coulomb potential when the electron and the hole are 
in the ground state, i.e. ip — > V'o.Oi (V 0,0 ) and in both the 
first and second excited state, i.e. ip — > Vo, l (V 0,1 ) and 
if) — > Vo,2 (l^ ' 2 ) respectively (the single particle densities 
are shown in Fig. 121) • Notice how all effective potentials 
in Fig. ^ show the —l/z behaviour for z — > oo and that 
the effective potential is finite for z = 0. 



V eS (z = 0) = -2.5961, 
d% s (z) 



dz 

1 d 2 V eS (z) 
2 



4.1967, 



dz 2 



-5.6140, 



(14a) 
(14b) 

(14c) 



2 = 



for the single particle ground state. In a similar way, we 
get a = -2.1126 (a = -1.9774), b = 3.1031 (b = 3.0608) 
and c = -4.5503 (c = -5.4277) for the first (second) 
excited state. This polynomial expression, Eq. (|13fl . is 
within 1 % of the numerical result for z/R < 0.14. 

Next we will look for an analytical expression that 
fits the whole range of z- values. Therefore, we have 
approximated the full effective potential by Pade ap- 
proximants. The first Pade approximation is given 
by 



^effW 



Po(z) = v 
Pi(z) w+\z\'- 



(15) 



where Pi(z) is a polynomial of z— th order. If we 
take into account that V^gf(0) = v/w and remember 
that for z — > oo the effective potential behaves as a 
normal Coulomb potential (—1/z), we can conclude that 
v = — 1 and calculate that w = 0.3852 (ground state), 
w = 0.4734 (first excited state) and w = 0.5057 (second 
excited state). 

This Pade approximation and the result of the cal- 
culations are compared in Fig. Ela) for the ground state 
and Fig. |3fb) for the first and second excited state. 
To improve this result we consider the next Pade approx- 
imation: 



B. Approximate analytical expressions 

Since the calculation of the effective Coulomb potential 
requires a considerable amount of computational time, it 
is highly desirable to have an analytical expression for it. 
Let us first consider the small z behaviour of the effective 
potential: 



bz 



cz 



(13) 



l\z\+8 



P 2 {z) z 2 +r 1 \z\+f3' 



(16) 



Again, there are certain conditions which impose 
restrictions to the values of the parameters. These 
conditions are V e ff(z) z i^f —1/z, which implies 
that 7 = — 1 and 8/(3 = V c q(z = 0), which leaves 
us with two fitting parameters. After fitting the 
data to the expression in Eq. (|16fl . we find rj = 
1.1288±0.0029, /3 = 0.4705±0.0015, S = -1.2215±0.0039 
for the electron and hole in the ground state, 
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FIG. 4: Comparison between the numerical effective electron- 
hole potential and the approximate one using the wave func- 
tion of a harmonic oscillator for three different sets of quan- 
tum numbers of electron and hole. 



r/ = 1.728 ± 0.017, [3 = 0.899 ± 0.011, S = 1.900 ± 0.023 
for the electron and hole in the first excited state and 
77 = 2.749 ± 0.061, (3 = 1.559 ± 0.039, S = -1.9774/3 = 
3.083 ± 0.077 for the electron and hole in the second 
excited state. The fitted second Pade approximation 
is also shown in Fig. [3] and we obtain an excellent fit 
(i.e. within 1.5%) for all states. 

If the confinement potential for the electron and 
hole is a parabolic potential, i.e. V(x, y) — mui 2 (x 2 +y 2 ), 
it is possible to perform all integrals in Eq. Ijllfl analyt- 
ically. The analytical result for the effective potential 
was obtained in Ref. [la and reads 



FIG. 5: Binding energy as a function of the parameter a = 
a* B /2R. The solid curves are calculated by using the effective 
potential V ' (electron and hole in the ground state). The 
symbols represent the fit which is given by Eq. 1191 . Inset: 
The same, but now for the effective potential V ' (electron 
and hole in the first excited state). 



using the finite difference technique. Fig. [S] shows the 
energy of the Coulomb interaction of the three lowest 
exciton energy levels as a function of the dimensionlcss 
parameter a — a* B /2R for two cases of the effective po- 
tential. 

The binding energy can also be fitted to an analytical 
expression now. Values of a < 0.1 are not realistic, since 
this would correspond to a large radius and a strong con- 
finement approach would no longer be valid. Therefore, 
the following fit is suggested for a values exceeding 0.1, 



"*<*>-- ( 3 ) m 



l-erf 



( Igj 



\V2lo/R 



Jz| 2 /2(/o/-R) 2 



(17) 



with the oscillator length Z = (h/mujo) 1 ^ 2 . It is remark- 
able that this potential approaches our effective exciton 
potential very closely if we define Iq to be 



lo/R 



^cff(O) 



7T\V2 

2, 



(18) 



i.e. the expression for the oscillator length extracted from 
Eq. (jT7|) by putting z = 0. It is possible now to calcu- 
late numerically Iq/R for the different effective poten- 
tials: l /R = 0.4828 (ground state), l /R = 0.5933 (first 
excited state) and Iq/R = 0.6338 (second excited state). 
Using these results in Eq. I|17|). we get an expression for 
the effective exciton potential as function of the inter- 
particle distance. The approximation and the calculated 
data can be seen in Fig. 01 and show a good agreement 
for ground state and first excited state (fit within 0.7% 
of the data), while for the second excited state the ap- 
proximated result is still reasonable (fits within 3.2% of 
the data). 



IV. EXCITON ENERGY 



1 + ba T 



(19) 



For the ground state configuration of the elec- 
tron and hole, we find that f = -2.145 ± 0.020, 
b = -1.093 ± 0.022, t = 0.5385 ± 0.0040, and for the 
first (second) excited states of the exciton binding en- 
ergy we find £ = -1.287 ± 0.031, b = 5.42 ± 0.16 
and r = 0.945 ± 0.013 (£ = -1.211 ± 0.016, 
b = 10.23 ± 0.14 and r = 0.8830 ± 0.0034) (see 
Fig. [5] for the fits). Similarly, for the configuration 
where electron and hole are in the first excited state, we 
find that £ = -1.714 ± 0.015, b = -0.967 ± 0.020 and 
r = 0.5449 ±0.0042 (ground state), £ = -1.0525 ± 0.028, 
b = 4.56±0.16, r = 0.940±0.016 (first excited state) and 
£ = -0.965 ± 0.016, b = 8.47 ± 0.15, r = 0.8944 ± 0.0051 
(second excited state). 

From Eq. (|12fl it is clear that for small a— values, 
there is a relatively smaller contribution of the kinetic 
energy term ~S7 2 ip- As a increases, the energy of the 
particle will also increase, because of an increase in 
the contribution from the kinetic energy term. For 
large values of a the energy levels become closely 
packed to the Ec = level. In fact, because of the 
1/z behaviour of the exciton potential for lar ge z , the 
integral / V(z) 1 / 2 dz = 00 and according to Ref. Ha there 
will be an infinite number of bound states. 



A. Numerical calculations 

When a single exciton is created by excitation, its energy 
is generally given by E — E g + E e + Eh + Ec, where 
Eg denotes the energy band gap, i? e (/i) is the single elec- 
tron (hole) energy and Ec is the (negative) energy of the 
Coulomb interaction between electron and hole. For a 
given material, we have to calculate now Ec, all other 
contributions to the exciton energy are known. There- 
fore, the ID Schrodinger equation [Eq. (|12fl ] was solved 



B. Comparison with analytical results 

In this section we present an analytical solution of the 
ID Schrodinger equation Eq. (|12|) . Similar calculations 
were done earlier by Loudon^, who solved the problem 
of the ID hydrogen atom (i.e. the ideal limit of the ID 
electron-hole system: infinitesimal wire cross-section and 
hard wall confinement) analytically. To avoid the diver- 
gence of the original potential at z = Loudon used 
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FIG. 6: Ground state (lower curves) and first excited state 
(upper curves) energy of the exciton when 1) using the first 
Pade approximation for the electron-hole interaction potential 
and 2) the numerical expression for the electron and hole in 
the ground state. The curve with symbols is the asymptotic 
behaviour, i.e. Eq. 12311 . 



different models, e.g. hard wall conditions at z = 0, 
i.e. cj)(z = 0) = 0. The ID Schrodinger equation, Eq. (JT3J, 
has in general to be solved numerically. We used the fi- 
nite difference technique. In the special case in which 
the effective potential is approximated by the first Pade 
approximant, Eq. (|15|) . it is possible to obtain an analy- 
tical expression for the wave function. After introducing 



the new variable y = 2y — Ec/a(z +w), we can rewrite 
Eq. 112(1 as the Whittaker equation: 



\n(w I ma) 



(odd states) 
(even states) 



(24) 



for m = 1, 2, 3, These correspond to the bound states 

of the exciton. This asymptotic result is shown in Fig. [B] 
by the square symbols and is clearly only valid for large a. 
Note that in the limit of a Coulomb interaction potential 
the excited states are twofold degenerate. The lowest 
bound state satisfies the special relation 



In 



0. 



w \ 1 
v a) 2v 

This approximation is also shown in Fig. 

V. MAGNETIC FIELD DEPENDENCE 



(25) 



0. 



(20) 



2yaJ -E c /a 



The wave functions are then given in terms of the Whit- 
taker function 



A. Single-particle properties 

When a magnetic field, directed along the wire, is ap- 
plied, the single-particle Hamiltonian for an electron 
(q = — |e|) or hole (q = |e|) can be written as 



0(y) = ^e^/ 2 (7(f- K ;2;y), 
is a normalisation constant and n 



(21) 



where A 

— v /2a\J — Ec/ct- The energy is then obtained by solving 
the following equations: 



My) 



dy 



(even states), 
cj)(y = 2w\]-Ec/a) = (odd states). (22) 



y=2wy/-E c /a 



We now obtain the energy of the Coulomb interaction 
semi- analytically (by solving the transcendent equation) 
and we compare it in Fig. with the numerically calcu- 
lated energy. Notice that for large a both energies are 
close to each other. But for small a-values (which cor- 
responds to large R) there is a underestimation of the 
exciton energy up to 15 %. This is not surprising, since 
the first Pade approximation leads to a rather poor fit to 
the effective potential (see Fig. 01a)). Nevertheless ana- 
lytical calculations are useful because they allow for an 
estimation of the energy without having to solve the dif- 
ferential equation (|12|l numerically. Furthermore it has 
to be noticed that, in the case of large a (which cor- 
responds to small radii R), it is now possible to find a 
simple analytical expression^ for Ec, i.e. 



H = -!- ( -iHV -qA) + V(r) 
2m 



(26) 



where V{r) is again the confinement potential of the wire. 
Due to the symmetry of the problem we may benefit from 
the symmetric gauge A — ^Bpeg for the vector potential. 
For the same reason, the wave function ip(p, 0) is separa- 
ble, and can be written as '^|}{p^ &) — R(p)e~ lW , where the 
z-dependence has been omitted for the sake of simplicity. 
The Schrodinger equation now reads 



dp 2 p dp p 2 



muj c 



, . /muj r \' z 2mE\ „ 

-hif) f + — ^°(27) 



where tu c — \e\ B /m is the cyclotron frequency. Introduc- 
ing the following parameters, k* = — sgn(q)lmLu c /h + 
2mE/h 2 and 1/1% = mui c /2h, and by putting R{p) — 
pMexp(-p 2 /2l 2 B )f(p), one gets 



tf"+((\l\ + l)-t)f- -(|Z| + 1) 



l 2 B k* 



/ = 0, (28) 



where the new variable t is given by t = p 2 /l 2 B . The 
solution is the confluent hypergeometric function! 18 . 



E c = -1 /4c 



(23) VVG°>0) = C n je- lW p 



where 



il8 n \ l \p 21 



» 1F1 



-a n x, \l\ + 1; 12 ) ' 



(29) 
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FIG. 7: Effective exciton potential for different values of the 
magnetic field (full lines) in which the electron and hole are 
in the ground state. The approximation (open symbols) of 
Eq. I|17|l can be used again to fit the calculated effective po- 
tential. For comparison, the effective potential for the system 
without magnetic field was also drawn (Ib/R = oo, dotted 
line). Note that for this figure, the electron and hole were in 
the ground state. 

where C n ^ is a normalisation constant and 

-a n ,l = \{\l\ + l)- 1 ^, (30) 

is a constant which is obtained by using the hard wall 
boundary conditions, i.e. ijj n j(R,0) = 0. The single par- 
ticle energy then becomes 

/ l + |Z|+sgn( g )^ 
hj n .i = Tuaj c I a„; H - I . (31) 



FIG. 8: Effective exciton potential at z = as function of the 
magnetic length Ib (the magnetic field). The fit is to Eq. (I35P 
and Eq. ff7J 

interesting to determine the relation between Ib and 
t4ff(z = 0), since the value of V(z — 0) will enable us 
to calculate Iq/R (Eq. H18|) ~) and therefore to reconstruct 
the effective potential for the whole range of z— values 
for an arbitrary magnetic field. Fig. [S] shows this curve 
for both electron and hole in the ground state and in 
the first excited state. 

This figure can be understood as follows. Large 
values of Ib correspond to the absence of a magnetic 
field and therefore the curve converges to the previously 
calculated value of V e Q t x (z = 0) for B = 0. For 
small values of Ib, Eq. I|17|l matches the calculated 
curve, because for high magnetic fields, the confinement 
becomes effectively parabolic. Replacing the oscillator 
length lo by Ib, the result for large magnetic fields is 
recovered: 



B. Effective exciton potential 



The total Hamiltonian is now given by 



Hint — — - 



2m 
1 



Tip \ 



Px e ,y c + \e\A\ + V e (x e ,y e ) 



2m h 
1 



Px h ,y h -\e\A) +V h (x h ,y h ) 

h 2 _ 



47re y/{x e - Xh) 2 + (Ve ~ Vh) 2 + Z 2 



(32) 



Again, we can reduce this equation to a ID effective 
Schrodinger equation similarly as for the B — case by: 
i) separating the motion of the two particles into the cen- 
ter of mass relative coordinates, ii) by using the adiabatic 
approximation for the exciton wave function, and iii) by 
introducing the same dimensionless units as in Sect. lIII Al 
This leads us to the equivalent of Eq. (|12f) . where the sin- 
gle particle energies are now given by Eq. (|31|l and the 
averaging of the effective potential has to be performed 
with the Kummer functions [Eq. (|29[) ] (instead of Bessel 
functions). Note that the confluent hypergeometric fun- 
tion \F\ is equivalent to the Kummer function M: 



L-Fi(a; b; z) — M(a, b, z). 



(33) 



Fig. [7] shows the effective potential for different values 
of the magnetic field (which corresponds to different 
values of the magnetic length Ib — \j2h/mu) c ). These 
potentials can be approximated again by the expression 
in Eq. (|17|) (see fit in Fig. 0). Therefore it is more 



eS,i B - 



^ = ^ = -\l-2Z 



/R 



(34) 



Putting together these results into a single formula, we 
obtain a fit for the effective potential as function of the 
magnetic field: 



1/0,0/ ,s ,yoo 

V eS ( z = °) = V eS, h 



(z = 0)- \l~J-e- ph 
W 2 1/R 



(35) 



where I /I 2 = l/l 2 + l/l B . Two parameters have 
been determined by fitting the data and are 
l /R = 0.0624 ± 0.0021 and pR = 8.570 ± 0.090 
for the electron and the hole in the ground state config- 
uration. Notice that basically we used the expression for 
the parabolic confinement Eq. 1|17|) an d by rescaling Iq 
to I, we include the effect of the magnetic field. 

By making a similar calculation to the one in Ref. IT3I 
but for the electron and hole in the first excited state, 
we find that 



eff,z f 



Xz = 0) 



11 
16 



1 



2 l B /R 



(36) 



Therefore we fitted the data for \n e ,l e ) = \rih,lh) = 
I to, I) = |0, 1) to a slightly different function, 



V^(z = 0) = V r [ 



r/0,l 
eff,i £ 



• v ' 16 V 2 Z/i? 



(37) 

where h/R = 0.0462 ± 0.0026 and pR = 10.73 ± 0.17. 
The numerical results together with the results of the fit 
are shown in Fig. [S] and agree very well. 
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FIG. 9: Ground state exciton binding energy for the electron 
and hole in the ground state as function of the magnetic length 
Ib = \/2h/eB for different values of a. Inset: the same but 
now for the ground state, first and second excited state with 
a = 1. 



C. Exciton energy 

Previously, we calculated the exciton energy as function 
of the dimensionless parameter a = a* B /2R, which is de- 
termined by the wire radius and the material parame- 
ters. In this case, there is another variable which can 
be changed and will have an effect on the energy: the 
magnetic length Ib- Considering the previous results of 
the system without magnetic field and knowing that the 
shape of the potential does not fundamentally change by 
applying a magnetic field (the effective potential well only 
gets deeper) , we know that also the energy as function of 
a will not give a fundamentally different result as before. 
Therefore we calculated the binding energy as function 
of the magnetic field (or the magnetic length Ib) for a 
few fixed values of a (a = 10, a — 1 and a — 0.1) and 
the result of these calculations can be seen in Fig. EI The 
main figure shows only the ground state binding energy, 
while the inset shows also the energy of the first and sec- 
ond excited state. The inset shows that the magnetic 
field influences the ground state more than the first and 
second excited states, which are almost constant. For 
high magnetic fields, which corresponds to Ib/R — > 0, 
we find that the effective potential of Eq. (fP7|l becomes a 
l/z-potentia! 15 . This means that we can use Eqs. 12311 - 
(|25J) again to determine the asymptotic behaviour of Ec 
as function of Ib- We find that the ground state energy 
Ec — > oo for Ib — > for all values of a, whereas for a = 1 
(inset of Fig.|5J the first and second excited state become, 
in the limit B — > oo degenerate with energy —1/4. 



VI. CONCLUSION 

In summary we have calculated the effective interaction 
potential and the binding energy for an exciton in a 
nanowire. The effective interaction potential was calcu- 
lated for different states of the electron and hole and we 
were able to present 'exact' numerical results. In order 
to reduce the amount of computational time for future 
calculations, these results were used to propose various 
analytical approximate expressions. We fitted the results 
to three different functions. A first Pade approximation 
gives the correct qualitative behaviour, but the fit was 
not optimal. Nevertheless this expression was useful to 
perform further analytical calculations to estimate the 
binding energy of the exciton, especially for small val- 
ues of the wire radius R. A second Pade approximation 
as well as the result from the parabolic confinement re- 
sulted in very good fits. Furthermore the exciton binding 
energy was calculated both semi-analytically and numer- 



ically and analytical approximations were given. We also 
investigated the influence of a magnetic field and deter- 
mined an analytical expression for the effective potential 
for arbitrary magnetic field values. Finally we calculated 
the binding energy as a function of the magnetic field, 
where we found that the ground state binding energy is 
most affected by a variation of the magnetic field. 
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APPENDIX A: NUMERICAL CALCULATION OF 
THE EFFECTIVE POTENTIAL 

To obtain the effective exciton potential of Eq. (|11|) , we 
have to evaluate four-fold integrals. They can be rewrit- 
ten in polar coordinates as follows, using the same di- 
mensionless units as before: 



V eS( z ) = -R 4 ffff PePhdp e dp h dQ e d9 h 

x \^ nl (Rp c ,e e )\ 2 \4, n>l ,(Rp h ,e h )\ 2 



(Al) 



Due to cylindrical symmetry, we can perform the integral 
over the angle Oh analytically, yielding 



d6 h 



d9 e / 



2tt 



V / {Pe+ Ph) 2 + Z 2 \{Pe+Ph) 2 +Z 



^PePh 



(A2) 



where K(m) is the elliptic integral of the first kind. Using 
a polynomial approximation^ for the elliptic integral, 



K(m) 



4 

E 

i=0 



4 1 

a,im\ + y biTn\ In — 
t=o 1 
= A(mi) +B(mi)ln(l/mi), 

with mi = 1 — to, we find that 

V cS (z) = -2irR 4 f dp ePe \ij nl {Rp e ,6 e )\ 2 
Jo 



(A3) 



^ dphPh \tpn'l'(Rph,0h)\" 
\/{pe + Ph) 2 + Z 2 



A(mi) +B(TOi)ln- 



In a next step, we calculate the integral I\ over using 
numerical techniques. To obtain good convergence, we 
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use the logarithmically weighted method and consider and the remaining problem is the calculation of the coef- 
the following integral: ficients Ai and C^. These coefficients are given by 



I(Pe,z) 



/ dp h F(p h ,p e ,z) In 
Jo 



(Pe - Ph) 2 + z"' 



(Pe + Ph) 2 + 



We use the following transformation 



(A4) 



I( Pe ,z) = £ I dxF( 

.•_n Ja 



x + hi, p e , z) 



x In 



(x — (p e — hi)) 2 + z 2 
(x + (fie + hi)) 2 + z 2 



, (A5) 



where h denotes the step size and s the total number of 
steps. If we replace F(x + hi) by Fi + (-Fj+i + F{)x/h, we 
can write l|A5[) as 



A t \ f A 1 f (x-(Pe-hl)) 2 + Z 2 

MP - Z) = Jo dxh \ (x + (p e + M)) 2 + z 2 
r< \ 1 t a i ( (x-(Pe-hi)) 2 + z 2 

aiPe ' Z) = hi dXXln {(x + (p e + M)) 2 + z 2 



(A7) 



These integrals can be performed analytically which leads 
to the following results: 



Ai(p e ,z) = a(h - (p e - hi)) - a(h + (p e + h{)), (A8a) 
Ci(p e ,z) = h^ 1 [c(h - (fie - hi)) - c(h + (p e + hi))] 
-2r e , (A8b) 



I(p e ,z)=Y, F *MPe,z) + (F l+1 -F t )C l (p e ,z), (A6) with 



i=0 



a(y) = H i (y) + H i (h-y) + 2 

mv) 
c{y) 
Gi(y) 



arctan 



(I) 



arctan 



h-y 



y\n(y 2 + z 2 ), 

\ (y(2h - y) + z 2 ) Gi(y) + \{(h~ yf - z 2 ) d(h - y)2z(h - y) 



arctan 



(?) 



arctan 



h-y 



Mir + * )• 



Now we can write, using the logarithmically weighted where 
method, 



(A9) 



h = 



idphph \i> n 'v{Rph>Qh)Y 



\/(Pe + Ph) 2 + Z 2 



(A(mi) +B(mi)ln(mi)) 



Fi = F(hi) 



4hi\*l> n ., l .(Rhi,0 h )\' 



s-l r h 



(All) 



y/(p e + hi) 2 + Z 2 

^2 / dxF(x + hi) (A( mj + B(mi) ln(mi)) 
i=o "'o 

, The final step is to calculate then the integral over p e . 

(FihA( mi ) + F (At + Q-! - d) B(mi)) Similarly, we can also perform this last integration over 



+F s -K(m x ) + F s C 5 _iB(mi), 



p e . Using the above mentioned result of Eq. I|A10|) . we 
obtain the final expression for the effective potential in 
(A10) the z-direction: 

I 



^eff(^) 



2tt 

'n 2 



R 4 



s-l 



hJ2 h j\Tpm(Rhj,8 e )\ 2 [[Si(hj)h + T(hj)^A( mi )] - [S 2 (hj) +T(hj)C s - 1 M(m 1 : 

L 3=1 



+ ^hs\4> n i(Rsh,9 e )\ 2 [[5i(/»s)ft + T(/is)^A(mi)] - [5 2 (H + T(HC s -iB(™i)] 
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with 



Si(x) 
T{x) 
S 2 (x) 



E 1 4hi\ip n >i>(Rhi,9 h )\ 2 
77 , -N9 7711 
i=1 ^(ar + hi) 2 + z 2 

4s/i|^ n /,/(A»fe,e h )| 2 



+ sft,) 2 + 2 2 



4/ti 1^/(^/11, g h ) |" 
V(x + hi) 2 + z 2 



E 



(A12) 

(A13) 

■ Ci_! - a)B(mi). 

(A14) 
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